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ABSTRACT 

Beyond several AU, interactions among shocks and streams give 
rise to "merged interaction regions" in which the magnetic field 
is turbulent. The integral intensity of > 75 MeV/Nuc cosmic rays 
at Voyager is generally observed to decrease when a "merged 
interaction region" moves past the spacecraft and to increase 
during the passage of a rarefaction region. When the separation 
between interaction regions is relatively large, the cosmic ray 
intensity tends to increase on a scale of a few months. This was 
the case at Voyager 1 from July 1, 1983 to May 1, 1984, when the 
spacecraft moved from 16.7 to 19.6 AU. Changes in cosmic ray 
intensity were related to the magnetic field strength in a simple 
way. It is estimated that the diffusion coefficient in merged 
interaction regions at this distance is v'0.6 x 10 22 cm 2 /s. 

1 • Introduction . Variations in the intensity of galactic cosmic rays > 
75 MeV/nucleon near 11 AU on scales ranging from approximately a day to a 
year were found to be related to the interplanetary magnetic field in 
observations from Voyagers 1 and 2 from June, 1982 to August, 1983 
(Burlaga et al., 1985a). The long-term variation of cosmic ray intensity 
was related to the strength and separation of interaction regions. It 
decreased or remained relatively low when interaction regions were strong 
and closely spaced, and it increased when the interaction regions were 
weaker and widely spaced. These results are consistent with the idea 
that modulation is caused by diffusion in turbulent magnetic fields 
(Morrison, 1956) and the observation that modulation effects propagate 
outward from the sun at the solar wind speed (McDonald et al., 1981). A 
non-steady model of the 11 -year variation which incorporates diffusion 
and propagation of shells of disturbances has been constructed by Perko 
and Fisk (1983). 

The purpose of this paper is to extend the analysis of Burlaga et 
al . (1985a) to examine the relation between cosmic ray intensity and the 
interplanetary magnetic field from July 1, 1983 to May 1, 1984, using 
Voyager 1 data. During this interval the spacecraft moved from a 
heliocentric distance of 16.7 AU to 19.6 AU and from a heliographic 
latitude of 20° to 23°. 
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2. Observations . The intensity of cosmic rays > 75 MeV/nucleon, measured 
by the CalTedh/University of New Hampshire/Goddard Space Flight Center 
experiment is shown in Figure 1, together with the strength of the 
interplanetary magnetic field measured by the GSFC magnetometer (N. Ness, 
Principal Investigator). The magnetic field strength was^normalized by 
the spiral magnetic field strength, B = 4.75 x (1 + r ) /r . Thus, 
the magnetic field fluctuations in Figure 1 are perturbations on the 
large-scale average magnetic field, which are produced by dynamical 
interplanetary processes. 

Figure ) shows a 
correlation between 
B/B and the cosmic 
ray p intensity C: C 

generally decreases 
during the passage of 
an interaction region 
or merged interaction 
region (B/B > 1) and 
it increases during 
the passage of a rare- 
faction region (B/B < 

1). A similar de- 
lation was observed 
near 11 AU by Burlaga 
et al. The net effect 
in the period shown in 
Figure 1 is an overall 
increase in cosmic ray 
intensity. The cosmic 
ray intensity profile 
is the result of a 
competition between 
the effects of 
interaction regions and those of rarefaction regions. In this interval, 
the interaction regions were relatively weak (B/B ) m < 2,5 and the 
corresponding decreases in C were relatively small, while the rarefaction 
regions were large in extent, allowing more than enough time to recover 
from the individual decreases. 

3. Model . Burlaga et al. (1985a) modeled the variation of cosmic ray 
intensity observed near 11 AU from July, 1982 to August, 1983 with the 
following set of equations: 


0 65r 


VOYAGER 1 


8 

z 0 55 

3 
O 
o 


INTEGRAL RATE 
NUCLEI >75 MeV/NUC 
24 HR AVG 


“l 1 r 


4 / 


045 h 


3 3 1 - 


4 







24 HR AVG 


1 j_ 

AUG 1 OCT 1 DEC 1 FEB 1 APR 1 


OCT 1 
1983 


FEB 1 APR I 

1984 


Fig. 1 


Intensity of cosmic rays > 75 
MeV/nuclei versus time (top) and magnetic 
field strength relative to the spiral 
field strength versus time (bottom). 


dC 

dt 

= - D (B/B -1 ) 
P 

when B/B > 1 
P 

(D 

dC 

= R 

when B/B < 1 

(2) 

dt 


P - 


are 

constants and B/B 

is the measured magnetic 

field 


where D and R 

strength as a function of time. UsSng the 24-hour averages of B/B (t) 
shown in Figure 1, choosing D = 0.004 (counts/sec/day) and R = 0.002 
(counts/sec/day) and taking the initial value of C at July 1, 1 983* these 
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equations were inte- 
grated to obtain the 
model cosmic ray in- 
tensity profile shown 
in Figure 2. The ob- 
served 24-hour aver- 
ages of C(t) for > 75 
MeV/nucleon cosmic 
rays are also plotted 
in Figure 2, for com- 
parison. Note that 
the model curve has 
been plotted with a 
constant offset of 
0.03 (counts/sec/day) 
for the sake of clari- 
ty, but its initial 
value can be chosen to 
be identical to the 
observed value for the 
integration of (1) and 

(2). Between March 10 and March 20 there was a decrease in cosmic ray 

B/B 6 "^ this 1 time 8 (see j£ edict ® d ’ b ® cause there was no large increased 
! ! gUre 1} ' Except for this anomaly, the model 

provides a very good approximation to the cosmic ray observations 
throughout the 10— month interval 



Fig. 2 The result of a model (see text) 
together with a plot of the observed 
counting rates. 


A theoretical basis for the model given by (1) and (2) has recentlv 
assumptions^ a " d U * <,985> ' f °“"< «>■* »„dar^ oSS.ET ' 


dC V , . V 

dt = _ L (C " C o ) T C oB 


SB 


(3) 


diffusion coefficient, V is the solar wind speed. L is a 
characteristic length and «B is a measure of the fluctuations of B. It 

ha ^ b !r\ 0bserved tb ?^ 3 6B/B “°*3 between 1 AU and 5 AU, varying slowly 

Jake 6B/B an - e o aS 3 Z* , 1982) ’ 30 U L reL 2 .STto 

( o\ hao (-R /B o in J3). In this case the second term on the RHS of 

> fcK ? same as (1), and one can use our value of D to estimatJ 
the diffusion coefficient K. With D = 0.004 (counts/sec/day), 

!^ 3 ’. C o = 0-55 counts/sec and «B/B = 0.3 one finds K = 0.6 x ioVs' 1 . 

we used SB/B n th ® Value chosen by Chih and Lee (because 

efficient- L A k* 3 d’ !* Ut lb is close to the value for the diffusion co- 
1 th by Perk ° and Flsk 0983 5 t0 describe the 11 -year varia- 

takJne cc-Hf th6 term ° n the RHS of (3) with R in (2), and 

implies L? 4 auJ 3 COUnts/sec ’ one finds that R = 0.002 (counts/sec/day) 

In Figure 3 we show spectra of the components of B (upper curve) the 

ffH^JomnJLn f l0WSr ° U [ Ve) and the magnetic helicity times frequency 
Ma/i Qf l r” f r e f Tt aVerage Voyager 1 data from July 1, 1983 to 
freeH«m 198 ^v, i ! bhe fast Fourler transform method with 26 degrees of 
eedom, without detrending or filtering the data (Matthaeus and 
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Goldstein, 1982). Positive values of magnetic helicity are denoted as 
circles, negative values as triangles. Assuming that plasma is convected 
past the spacecraft at the mean solar wind speed in the interval, V» the 
frequency f corresponding to the correlation length L, is f Q = V/L, and 
this is shoSn by the arrow in Figure 3. 


The spectrum of 
power in the com- 
ponents pf-jj has the 
form f ^ expected 
for homogeneous turbu- 
lence in the range 4 * 

10"° Hz to 2 x 10~ D 
Hz, corresponding to 
periods from 14 hours 
to 3 days. This is 
consistent with 

earlier results show- 
ing that the fluctua- 
tions in interaction 
regions are turbulent. 

At periods between 4 
days and 15 days, the 
spectrum of power in 
the components of § 

was , which is 3 power spectra of the magnetic field (see 

probably either a text) . 

remnant of the spec- 
trum of fluctuations introduced at the source (Goldstein et al., 1984; 
Burlaga et al., 1985b) or evidence of an inverse cascade of magnetic 
helicity expected in fully developed MHD turbulence (Frisch et al., 1975; 
Montgomery and Matthaeus, 1981). If the former interpretation is cor- 
rect, its presence is another indication that turbulence did not develop 
sufficiently to modify the initial spectrum in this quasi-stationary flow 
system, despite the long time available for evolution, viz, ^18 AU/400 
km/s -► 78 days. 
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